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The observation of the F(4140) structure in — > J/ip (fyK^ decays produced in pp collisions 
at \/s = 1.96 TeV is reported with a statistical significance greater than 5 standard deviations. A 

fit to the J/xji (f> mass spectrum is performed assuming the presence of a Broit-Wignor resonance. 
The fit yields a signal of 19 ± 6(stat) ± 3(syst) resonance events, and resonance mass and width 
of 4143.4l^:g(stat) ± 0.6(syst) MeV/c^ and 15.3t6°i''(stat) ± 2.5(syst) MeV/c^ respectively. The 
parameters of this resonance-like structure are consistent with values reported from an earlier CDF 
analysis. 

PACS numbers: 14.40.Gx, 13.25.Gv, 12.39.Mk 



The existence of exotic mesons beyond qq has been dis- 
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cussed for many years [l|, but evidence for such mesons 
has not been clearly established. The recent discov- 
eries of states with charmonium-like decay modes [3- 
0] that do not fit into the overall charmonium system 
have introduced challenges to the conventional qq meson 
model. The possible interpretations beyond qq such as 
hybrid (qqg) and four-quark states {qqqq) have revital- 
ized interest in exotic mesons in the charm sector ■ 
Recently, the CDF collaboration has reported evi- 
dence for a narrow structure near the J /ip </> threshold 
in i?+ — )• J /ip (j)K~^ decays produced in pp collisions 
at y/s ~ 1.96 TeV 12|. Charge conjugation is implied 
throughout this letter. Since the mass of this state, 
termed y(4140), is well beyond the threshold of open 
charm-pair production, the expected branching fraction 
into this channel for conventional charmonium is tiny. 
The structure is the first observed charmonium-like ob- 
ject decaying into a pair of quarkonium states (cc and ss) 
with a relative narrow width, a possible signature for an 
exotic meson 0, [III [3, Q ■ The Belle collaboration has 
searched for this J/tl^cj) structure without a firm conclu- 
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sion near the J/ipfj) threshold 

In this Letter, we report a further study of the struc- 
tures in the J/ipcj) system produced in exclusive ^■ 
J/il!(j>K~^ decays with J/4> — 4> ~^ K^K~ 
reported in Ref. This analysis is based on a sam- 

ple of pp collision data collected by the CDF II detector 
with an integrated luminosity of 6.0 fb~^. This analysis 
includes the data used in, and supersedes the results of 
Ref. 0. 

The CDF II detector has been described in detail else- 
where [l6| . The important components for this analysis 
include the tracking, muon, and time-of-flight (TOF) sys- 
tems. The tracking system is composed of a silicon-strip 
vertex detector surrounded by an open-cell drift chamber 
system (COT) located inside a solenoid with a 1.4 T mag- 
netic field. The COT and silicon-strip vertex detector are 
used for the measurement of charged-particle trajectories 
and decay locations. In addition, the COT provides ion- 
ization energy loss information, dE/dx, used for kaon 
discrimination, while the TOF system provides comple- 
mentary kaon discrimination information. The central 
muon identification system is located radially outside the 
electromagnetic and hadronic calorimeters and consists 
of two sets of drift chambers and scintillation counters. 
The central detector covers the pseudorapidity region 
\r]\ < 0.6 and detects muons with px > 1.4 GeV/c [nj, 
and the outer part covers the region 0.6 < \r]\ < 1.0 and 
detects muons with pT > 2.0 GeV/c. 

In this analysis, J/ip — >■ events are recorded us- 

ing a dedicated three-level dimuon trigger. The first trig- 
ger level requires two muon candidates with two COT 
tracks that extrapolate to track segments in the muon 
detectors. The second level applies additional kinematic 
requirements to the muon pair candidate. The third level 
requires the invariant mass of the pair to be within 



the mass range of 2.7 to 4.0 GeV/c^. The trigger require- 
ments are confirmed offline. 

We apply the same requirements described in the pre- 
vious analysis fj] to the current data. We form i3+ ^■ 
J/^(j)K~^ candidates by combining a J/4' ^ can- 
didate, a — >■ K^K^ candidate, and an additional 
charged track, which are consistent with originating from 
a common point. The three hadronic tracks must be 
identified as kaon candidates by using a log-likelihood 
ratio estimator. This quantity reflects how well a can- 
didate track can be positively identified as a kaon rel- 
ative to other hadrons with its dE/dx and TOF infor- 
mation and must exceed 0.2 [l^. The reconstructed 
masses of the J/ijj and (f) meson candidates must lie 
within 50 and 7 MeV/c^ of their nominal values, respec- 
tively. In the final B"*" reconstruction the mass is 
constrained to the known J/ ip mass [H, and the B+ can- 
didates must have > 4 GeV/c. In addition, we require 
Lxy{B+) > 500 //m for the B+ -J> J/il}4>K+ candidate, 
where Lxy{B'^) is the projection onto pt{B^) of the vec- 
tor connecting the primary interaction point, determined 
for each event using prompt tracks, to the reconstructed 
B^ decay point. 

The invariant mass spectrum of the selected J/ tp (j)K^ 
candidates is shown in Fig. [ija) . It is fit with a Gaussian 
signal function with its root-mean-square (RMS) width 
fixed to the mass resolution of 5.9 MeV/c^ obtained from 
Monte Carlo (MC) simulation fl^ and mean fixed to the 
nominal B~^ mass [l|] and a linear background function. 
The 5+ yield is 115 ± 12(stat) events, a 53% increase 
over the previous analysis. This increase in yield, for 
an integrated luminosity increased by a factor of 2.2, is 
reduced by a trigger rate limitation at the higher in- 
stantaneous luminosities for the later data-taking pe- 
riod. The yield increase in the complementary mode 
B+ J/tI>tt+tt-K+ is 51.8±2.4%, consistent with the 
yield increase in 5+ J/ip(j)K'^ channel. 

We then select B^ signal candidates with a mass 
within ±3 RMS (±17.7 MeV/c^) of the nominal B+ 
mass. Events with a mass within [—9, —6] RMS or 
[-1-6, -1-9] RMS of the nominal B^ mass are called B side- 
band events. They are normalized into the signal 
region assuming a linear background distribution. The 
J/ip signal, checked by removing its mass constraint, 
contains almost no background. Figure [Hb) shows the 
invariant mass distribution of the K~^K~ pairs from 
J /ip K~ candidates inside the B mass window 
and in the B sidebands before applying the restriction 
on the K^K~ mass. The clear (p signal inside the B 
mass window and almost featureless K^K~ mass dis- 
tribution in the B sideband indicate that the 5+ 
J /ip K'^K~K'^ final state is well described as J/ip<pK^ . 
In none of the candidate events do both K~^K~ combina- 
tions from the three-kaon final state fall into the (p mass 
window. 

Fig. [5] shows the mass difference AM = 
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FIG. 1: (a) The mass distribution of J/ip 4>K^; the solid blue 
line is a fit to the data with a Gaussian signal function and 
linear background function, (b) The K^K~ mass distribu- 
tions inside the B mass window (black solid) and in the B 
sidebands (red dotted). 



m,{fi~^li~K~^K~) — m{ii^fj,~) for events in the 
mass window. Events from reference [l^ and from new 
data are shown in (a) top and bottom. In the F(4140) 
signal region (AM < 1.07GeV/c^), the new data agree 
within la of the expectation (6 events compared to 7.3 
expected). Over the entire examined region the two 
data sets are consistent at the 7% probability level. We 
have investigated the consistency of particle ID for the 
two data sets using the i?+ — > J/ip K'^ channel and 
see no discrepant effects. In (b) and (c), we display 
AM distributions for the events in the B signal and 
sideband in the combined data sample. We restrict our 
study to events with AAf smaller than 1.56 GeV/c^ to 
avoid appreciable combinatorial backgrounds from 
misidentified B° V'(25')0 { J/ipTr+Tr~)(t> decays [l^. 
An enhancement is observed near the J/^l> 4> threshold 
from the signal while there are no events in the 
AM range below 1.1 GeV/c^ from the combinatorial 
background estimated from B sideband events. 

We model the observed threshold structure by an 5*- 

con- 



wave relativistic Breit-Wigner (BW) function [21 



voluted with a Gaussian resolution function with the 
RMS fixed to 1.7 MeV/c^ obtained from MC. Three-body 
phase space [H is used to describe the background shape. 
There is still a small i?° contribution (3.3±1.0 events) 
in the AM distribution up to 1.56 GeV. The MC shape 
of the Bl contribution is normalized to this area and 
added to the three-body phase space. The parameters 
from an unbinned likelihood fit to the AA/ distribution, 
as shown in Fig. [2lb) , are given in Table H To test the 
hypothesis that the structure has zero width (weak de- 
cay), we also fit the AM distribution to a zero- width 
peak, using a single Gaussian with RMS given by the 
expected mass resolution (1.7 MeV/c^), plus phase space 
background. The statistical significance for a non-zero 
width determined by the likelihood ratio between these 
two fits is 3.7fT, favoring a strong decay (non-zero width) 
rather than a weak decay for this structure. 



FIG. 2: (a) The mass difference, AM, between n+ n~ K+ K~ 
and jj.'^ fi~ , in the B"*" mass window. Top-data from Ref. p^ . 
bottom-new data, (b) A fit to the combined data assum- 
ing y(4140) only, (c) A fit to the combined data assuming 
two structures. This fit, including the second peak, lowers 
the 3-body phase space background under the first peak and 
increases its yield and significance with negligible effect on 
its resonance parameters. The shaded histogram is the data 
from the B sideband. The dotted blue curve is the predicted 
background contribution, the dash-dotted black curve is the 
predicted contamination, and the solid red curve is the 
total unbinned fit. 



The combinatorial background contains primarily 
misidentified 4> candidates, as can be seen in Fig. [1] (b). 
These two tracks with a (p-like mass will be combined 
with a real J/il', and an additional kaon candidate, all 
having a common vertex and forming a B mass. We 
model this component with phase space. To check this 
assumption, we performed several studies in which we re- 
laxed cuts that would not influence the mass-difference 
distribution of events from the B mass region: loosened 
vertex requirements or loosened Lxy cuts. These studies 
show that the combinatoric background from the B side- 
band region is consistent with 3-body phase space. We 
can now conclude that the flat background hypothesis 
used in the previous paper [l^ was overly conservative. 

We determine the signiflcance of the structure at the 
J/tp <j> threshold based on simulation. We generated 
8.4 X 10'' mass spectra (119 events for each, correspond- 
ing to the number of observed events) drawn from a 
three-body phase-space-like distribution, and search for 
the most significant fluctuation in each spectrum in the 
mass range of 1.02 to 1.56 GeV/c^, with widths in the 
range of resolution up to 120 MeV/c^ 1^. We evaluate 
2Aln£ = —2ln{£,o / Lmax) value for each generated spec- 
trum, where Co and Lmax are the likelihood values for 
the null hypothesis fit and signal hypothesis fit. Both 



fits use three-body phase space to describe the back- 
ground. There are 19 generated spectra with a 2Aln£ 
value greater than or equal to the value (34.9 obtained 
in the data assuming the y(4140) structure only |23l] ) 
obtained in the data. The resulting p-value, taken as 
the fraction of the generated spectra with a 2Aln£ value 
greater than or equal to the value obtained in the data, 
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is 2.3 X 10 ^, corresponding to a significance greater than 
5.0a [H. 

The mass of this enhancement is 4143. 4+3 g MeV/c^ 
after adding the J/ijj mass Q to the AM calculation. To 
study the systematic uncertainties of the mass, width, 
and yield, we repeat the fit to the AM distribution us- 
ing a non-relativistic BW and P-wave relativistic BW 
for signal. Other systematic uncertainties were also con- 
sidered, including assuming the existence of the second 
structure, varying the mass resolution and i?^ component 
amplitude, as well as the systematic uncertainty due to 
the particle identification. The resulting systematic un- 
certainties for the measured quantities are shown in Ta- 
ble H 



TABLE I: The fit results for J/tp(j} resonance near thresh- 
old. The first uncertainty is statistical, and the second one is 
systematic. 



AM [MeV/c^] 


Width [MeV/c^] 


Yield 


1046.7^3;; ±0.6 


15.3+6"!* ± 2.5 


19±6±3 



The relative trigger and reconstruction efficiency 
e{B+ Y{4U0)K+) x e(y(4140) J / (/)) / e{B+ 
J/%p(j>K'^) is determined to be 1.1, using an 5'- wave 
BW with mean and width values determined from data 
to represent the y(4140) structure and a three-body 
phase space kinematics for the J/-0 (f)K'^ de- 

cay. Thus the relative branching fraction Brei = 
B{B+ y(4140)/f+) X S(r(4140) ^ J /i; (j)) /B[B+ 
J/tl;<j)K^) including systematic uncertainties is 0.149 ± 
0.039(stat) ± 0.024(syst). 

An excess above the three-body phase space back- 
ground shape appears at approximately 1.18 GeV/c^ in 
Fig. [5] (b). Since the significance of y(4140) is greater 
than 5cr, we assume the existence of the F(4140) with 
the parameters given in Table U and background given 
by three-body phase space, and we test for the exis- 
tence of a possible structure around 1.18 GeV/ as shown 
in Fig. [5] (c). The signal PDF for the second struc- 
ture is an iS-wave relativistic BW function 21] convo- 
luted with a Gaussian resolution function with the RMS 
fixed to 3.0 MeV/c^ obtained from MC. For the second 
structure —2lii{Co / Cmax) is 16.8, where Cq and Cmax 
are the likelihood values for the null hypothesis fit as- 
suming the F(4140)-only and signal hypothesis fit as- 
suming the F(4140) and a second structure near AM 
~ 1.18 GeV/c^. The p- value determined by a simula- 
tion similar to the F(4140) investigation is 1.1 x 10^^, 
which corresponds to a significance of 3.1cr. The fit 
returns a yield of 22 ± 8 events, a AM of llTT.Ttg;^ 
MeV/c2, and a width of 32.3l?^;| MeV/c^ for the struc- 
ture near AM ~ 1.18 GeV/c^. Refitting the AM distribu- 
tion with a second structure produces negligible changes 
in the mass and width of the y(4140). The yield of 



the y(4140) increases by one event. We evaluated the 
systematic uncertainties for the second structure in the 
same way as for the y(4140) structure and found sys- 
tematic uncertainties of 1.9 MeV/c^ for the mass and 
7.6 MeV/c^ for the width. The mass of the second struc- 
ture is 4274.4t^:^(stat)±1.9(syst) MeV/c^ after including 
the world-average J/ip mass. 

In summary, the increased B^ — > J/ijj (j)K^ sample at 
CDF enables us to observe the F(4140) structure with 
a significance greater than 5cr. Assuming an 5'-wave rel- 
ativistic BW, the mass and width of this structure are 
measured to be 4143.4;|^(stat) ± 0.6(syst) MeV/c^ and 
15.3t6V(stat) ± 2.5(syst) MeV/c^, respectively. They 
are consistent with the previous report (m = 4143.0 ± 
2.9(stat) ± 1.2(s3^ MeV/c^, F = 11.7+|g(stat) ± 
3.7(syst) MeV/c^ [l4|. The relative branching fraction is 
determined to be Brei = 0.149±0.039(stat)±0.024(syst). 
Light meson vector-vector threshold enhancements have 
been seen [22[ . We do not know of any non-exotic mecha- 
nism for producing a threshold enhancement involving a 
pair of heavy quarkonium states, but we cannot exclude 
the possibility. We also find evidence for a second struc- 
ture with a mass of 4274.4t6;7(stat) ± 1.9(syst) MeV/c^, 
a width of 32.3ti5;3(stat) ± 7.6(syst) MeV/c^ and a yield 
of 22 ± 8 events. The significance of the second structure 
is estimated to be approximately 3.1cr. 
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